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1 Introduction 

Currently there is a great deal of evidence for the existence of dark matter in 
the universe: in our galaxy, in other galaxies ai^in-galactic clusters, as much 
as 90% of matter in the universe may be darffla&cl. There are many condi- 
dates for dark matter in astronomy such as jupiters, brown dwarfs, blackholcs, 
etc. However, not all of dark matter can be baryonic and thus a number of 
particle physics candidates have been suggested. These include massive neutri- 
nos, axions, supersymmetry candidates as well as exotic candidates such as 
topological defects. A convenient parametrizaton of matter in the universe is 
given by the quantity, fl=p/p c , where p is the matter density of the universe 
and p c = 3 H 2 /8ttGn is the critical matter needed to close the universe, and 
H is the Hubble parameter which is generally parametrized by H=100 h km/s 

Mpc. Experimentally the current evaluations of h lie in the interval 0.4 ~ h 
^0.8. An attractive hypothesis theoretically is that of an inflationary universe 
with Q = 1. The quantity that enters in theoretical analyses of dark matter 
is £lh 2 . The COBE data indicates that there are more than one component 
to the non-baryonic dark matter, e.g., a hot component comprised of particles 
which would be relativistic at the time of galaxy formation and a cold com- 
ponent comprised of particles which would be non-relativistic at the time of 
galaxy formation. The hot component could be a massive neutrino and the 
cold component could be either axions, or a susy candidate. Actually, in ad- 
dition to the cold-hot dark matter fit to the power spectrum data observed 
by COBE and other experiments there are other possibilities as well, such as 
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inclusion of a cosmological constant. We list below some examples: (i) CHDM 
(cold-hot DM with %o = 0.7, n HDM = 0.25, n B = 0.05 (baryonic DM), h = 
0.5) giving f^o h 2 = 0.18; (ii) ACDM (cosmological constant DM with Q^o = 
0.25-0.35, 0^= 0.05, Q A = 0.6-0.7, h =0.8) with fi^o h 2 = 0.16-0.23.; (hi) 
tCMD (tilted DM with fi^o = 0.95, fi B = 0.05, h = 0^5-0.6, n = 0.8) giving 
J7^o h 2 — 0.24-0.34. In the analysis of this review we shall assume that the 
parameter space is limited by 

0.1 < O^o h 2 < 0.4 (1) 

We will discuss later on the effect of varying the allowed range of f^o h 2 . 

Theoretically the analysis of f2^o h 2 is interesting as it involves several 
different fields of physics: particle physics enters via the fundamental inter- 
actions of the dark matter candidates, cosmology enters as the relic density 
involves both the gravitational constant Gn and the Hubble parameter H, and 
statistical physics enters via the appearance of the Botzmann constant k in the 
analyses of relic density in the early universe. Further, nuclear physics enters 
as one considers specific detection techniques such as scattering of dark matter 
particles from nuclei. 

In this review our focus will be to explore in detail the supersymmetric 
candidates for cold dark matter. In the Minimal Supersymmetric Standard 
Model (MSSM) there are two candidates which are neutral and have odd R 
parity: these are the neutralino and the sneutrino. If we assume that one of 
these is the lowest mass supersymmetric particle (LSP) and further assume 
that R parity is conserved then, one of the candidates could be a CDM candi- 
date. Unfortunately there are 110 arbitrary parameters in MSSM, and thus the 
masses of the supersymmetric particles are all arbitrary. Thus computations 
of supersymmetric dark matter in the framework of MSSM are by their very 
nature hi&hly unpredictive. In contrast, in the minimal supergravity grand 
unificationp'El one can describe low energy physics in terms of just four param- 
eters. Thus the theory is very predictive. 

A topic of great interest is the laboratory detection of dark matter. In 
this context the scattering of dark matter particles from nuclei appears very 
promising and we will devote considerable amount of discussion to this topic 
in this paper. The outline of this paper is as follows: In Sec 2 we discuss dark 
matter in supergravity unification; in Sec 3 a discussion of the relic density 
analysis is given; in Sec 4 we discuss constraints on dark matter analyses; Sec 
5 is devoted to a discussion of the detection of dark matter; in Sec 6 we give an 
analysis of event rates for the case of universal soft SUSY breaking parameters; 



2 



in Sec 7 we discuss the effects of non-universal breaking; in Sec 8 an analysis 
of event rates with non-universal soft SUSY breaking is given and in Sec 9 we 
give conclusions. 

2 Dark Matter in Supergravity Grand Unification 

As discussed above the MSSM is very unpredictive. The reason for this is 
that there is no phenomenologically acceptable mechanism for the breaking of 
supersymmetry in the MSSM. Consequently one must add by hand the soft 
SUSY breaking parameters and there are as many as 110 arbitrary parameters 
that one can add in the soft SUSY breaking sector. For this reason we shall 
work here within the framework of supergravity grand unification which has a 
much more limited parameter space. The models we have in mind are those 
where suspersymmetry is broken in the hidden sector and communicated to the 
physical sector by gravitational interactions, and we consider here the minimal 
model. We assume that the basic theory in the physical sector is a grand 
unified theory with a GUT group G which breaks at a scale Mqut — 10 16 to 
SU(3) x SU(2) x U(l). After the breaking of supersymmetry and the GUT 
symmetry, and after integration over the fields of the hidden sector and the 
superheavy fields, one finds the effective potential in the SUSY breaking sector 
to be of the form: 

V S b = m E 

a z a z a + {A W Y + B oi i H x H 2 + h.c) (2) 

Lmass = — m l/2-^ Q ^ Q (3) 

where {z a } are the scalar fields , A a the gaugino fields, n is the Higgs mixing 
parameter and Wy is the cubic Yukawa coupling which enters in the effective 
superpotential,W, of the low energy theory as follows: 

W = fioH^ + W Y + -J-W< 4 > (4) 

Here contains any quartic non-renormalizable couplings which may lead 

to proton decay. At this level the theory contains the four SUSY breaking 
parameters mo, toi/2, Aq, and Bq in addition to the Higgs mixing parameter 
/io, and the GUT parameters Mq and qq. 

The dynamics of supergravity grand r mpification breaks the electro- weak 
symmetry dynamically by radiative effectUutl Further, one typically uses the 
radiative breaking of the electro- weak symmetry to compute the Higgs mixing 
parameter fi by fixing the Z mass, and one eliminates the parameter B(which 
is the value of the parameter Bq at the electro-weak scale) in favor of 
tan/3 =< H 2 > / < Hi > such that 
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, it? — rintan 2 3 1 „ , 9 —2Bu 

2 = pi ^2 _ _ m|; sin2/3 = — ^- — ~ (5 

P ian 2 /3 - 1 2 z ' M 2/i 2 + ^i 2 + fi 2 w 

Here /U 2 = m 2 ^ . + Ej where mf^ . is the running Hi mass at scale Q w Mz and 
Si are loop corrections. The m^. are given by 

= m l +rn\ /2 g{t) (6) 



™ff 2 = ™i/2 e W + A m 1/2 f{t) + m 2 h(t) - k(t)A 2 Q (7) 

where the form factors e,f,g,h,k are defined in ref(10) and t=\n(M G / Q 2 ) . One 
may show that solutions exist to Eqs. (5), i.e. that SU(2)xU(l) is sponta- 
neously broken, if and only if at least one of the supersymmetry soft break- 
ing interactions are non-zero. Thus it is the supergravity interactions at Mq 
that give rise to the breaking of SU(2)xU(l) at the electroweak scale. The 
SU(3) x SU(2) x U(l) gaugino masses at the electro-weak scale are given by 



to; = (a i /a G )m 1/2 (8) 

In the above supergravity grand unification the sparticle mass spectrum is 
determined dynamically in terms of the four parameters, and thus the theory 
predicts what the LSP is. It is quite remarkable that over most of the parameter 
space one finds that the LSP is indeed a neutralino Xi- Thus supergravity 
predicts the x? to be the cold dark matter of the universe. The xi is of course 
an admixture of gauginos W3, B and Higgsinos H\, H2 and one may write 

Xi = niW 3 + n 2 B + n 3 H x + n 4 H 2 (9) 

where we use the convention that xi is the lightest neutralino and where 
are the components of the gaugino states and are to be determined by diago- 
nalizing the neutralino mass matrix. The neutralino mass matrix in the basis 
(W3, B, Hi,H 2 ) is given as follows: 



TO 2 
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b 
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TOl 


c 
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a 


C 





-A* 


b 


d 








(10) 



Here a = MzcosQwcosd, b — —Mzcos8wsin[3, c — —MzsinOwcosd and 
d = Mzsin9wsin3. In general one needs a numerical analysis to determine 
the value of the mass eigenvalues as well as the eigenvectors. In the limit 
/1 2 >> M§ , which holds in a significant part of the parameter space, it is poss- 
sible to obtain an analytic form for the indices n, as an expansion in 11/ Mz- 
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The results are given in appendix A. This expansion is useful in understanding 
the phenomena associated with the neutralino relic densities and event rates. 
Generally one finds that in this limit the neutralino is mostly a Bino with 
n,2 > 0.95 and often larger, while 713,714 and n\ are typically of first order, 
0{Mz I [i)- Because of this the Higgsino components are not totally negligible, 
and they are often quite significant e.g. 713 w 0.2. Hence one cannot neglect 
the Higgsino components even when the neutralino may be dominantly a Bino. 
This result has important implications in the analysis of event rates. 

Before proceeding further we wish to discuss the effect that a heavy top 
quark has on the electro-weak breaking. One of the major effects of a heavy 
top arises because it lies close to the Landau pole value which is given by 
771/ w 197sinf3. This occurs because the proximity of the top to the Landau 
pole drives the lighter stop towards a tachyonic limit. We exhibit this phe- 
nomenon in some detail. The mass of the lighter stop t\ is determined by the 
diagonalization of the following stop (mass) 2 matrix: 

m 2 R m t (A t - iicbnp) 
m t (A t — fJ.ctn/3) m| 

The lighter stop (mass) 2 eigen- value in the vicinity of the Landau pole is given 

m i = -llt + m U NP ) ^ 

where 

A R S A t - 0.6137775; D Q = 0.164[(mf ) 2 - (m t )} 2 /M 2 , (13) 

and where m 2 (NP) is a non-pole contribution. The condition that t\ not turn 
tachyonic limits severely the allowed range of A t . 

3 Relic Density Analysis 

An important quantity in the analysis of supersymmetric dark matter is the 
LSP relic density. After the time of the Big Bang, the primordial neutralinos 
begin to annihilate into ordinary matter, e-g., via the process Xi + Xi ~ * 
/ + /,..., where /, / are leptons, quarks, et&2. The annihilation is governed 
by the Boltzman equation : 

^ = -3Hn- <<rv> (n 2 - ng) (14) 
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where n = number/ vol is the Xi number density, no is the value of n at thermal 
equilibrium, H is the Hubble parameter at time t, a is the neutralino annihi- 
lation cross section, v is the relative neutralino velocity, and < av > stands 
for the thermal average. Assuming that the neutralinos are non-relativistic at 
the time of interest (a realistic approximation) one has 

/>oo />oo 

< av >= / dvv 2 crv exp[—v 2 /Ax]/ / dvv 2 exp[—v 2 / Ax\ (15) 
Jo Jo 

At large temperatures close to the Big Bang the Xi is m thermal equilib- 
rium with the background. However, the Xi decouples from the background 
when the annihilation rate becomes smaller than the expansion rate given by 
the Hubble term. This decoupling takes place at the "freezeout" temperature 
Tf. Integration of the Boltzman equation from the freezeout to the current 
temperatures gives the current relic density. One finds 



where 



n.oh 2 s* 2.48 x ltr 11 -±M — -r-^ (16) 

\Tj \2.73j J(x f ) { ' 



J (x f ) = / dx {av) (x)GeV- 2 (17) 



and Xf = kTf / . Now numerical analyses show that freezeout occurs for 
values of Xf ~ ^ which indeed corresponds to the neutralinos being non- 
relativistic at the time of the freezeout. The central quantity in the analysis of 
the relic density is J(xf). The computation of this quantity is quit subtle as it 
involves a double integral over (av) which in some regions is rapidly varying. 
An approximation that has generally been made is to expand av in a power 
expansion in v, i.e. av = a + b(v 2 / c 2 ) + • • •. With this approximation the 
computation of J(a;/) becomes trivial. However, it is knowro that such an 
expansion is not a good approximation when (av) is rapidly varying. This is 
precisely the situation that occurs in the ^Lcmihilation of the neutralinos when 
one is close to the Z and the Higgs poleoOij. In this case one must carry 
out a thermal averaging keeping the pole structure intact. To circumvent the 
problem of the double integration over the pole, one can carry out analyti- 
cally one of the integrations reducing J(x/) to a single integration which is 
computationally easy to evaluate. Details are given in Appendix B. 
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4 Constraints on Dark Matter Analyses 

The analysis of dark matter relic density is subject to many constraints. Since 
the analysis is carried out in the framework of supergravity grand unification 
one computes masses of 32 supersymmetric particles in terms of the basic 
SUSY breaking parameters (which are just four for the universal soft SUSY 
breaking case). The analysis of the mass spectrum is subject to following 
constraints: (i) color and charge conservation at the electro- weak scale, (ii) 
naturalness constraints, i.e., that mo,m§ < lTeV,(iii) LEP, CDF and DO limits 
on the mass spectrum. Further, the decay b — > s + 7 represents an important 
constraint on any new physics beyond the Standard—Model. Recently CLEO 
measured the branching ratio for B — » X s + 7 to bctll: 

BR(B -> X„i)) = (2.32 ± 0.5 ± 0.29 ± 0.32) x 10~ 4 (18) 

Combining the errors in quadrature gives BR(B— > Y s 7) = (2.32 ± 0.66) x 
10~ 4 .In the spectator approximation one can relate the above to the flavor 
changing process b — > s + 7 and one has 

BR{B^X sl ) _ r(6^s + 7 ) 



BR(B -> X c ev e ) T(b^c + e + v e ) 



R (19) 



where BR(B^ X c eP e ) = (10. 7± 0.5). The quantity R turns out to be a 
useful parameter as the factor (mjj) (which contains a significant experimental 
ambiguity) as well as some CKM factors cancel out in the this ratio. 

The b — ► s + 7 decay proceeds via loop corrections. For the case of the 
Standard Model the loop contribution involves a W-t exchange. For the su- 
persymmetric case one has in addition a,W — t, H~ — t and a Z — b exchange. 
We see then that the SUSY contributions enter on an equal footing with the 
Standard Model contribution. Thus b — > s + 7 is an important constraint on 
the SUSY parameter space. The supersymmetric contribution to the decay 
b — > s + 7 is governed by effective HamiltonianiS 

H eff = V tb V; s ^C 7 (M w )Q 7 (20) 

where F^ v is the electromagnetic field strength and Q7 = (e/2A-K 2 )mi ) SLC^ l 'biiF^ u . 
However, the computation of b — > s+7 must be carried out at the scale fi ~ rn;, 
where the decay occurs and thus one must use the renormalization group equa- 
tions to go from the scale \x « Myy to fi ~ mj. Using the above evolution, the 
value of R to leading order (LO) in QCD corrections is given by 

«=|^O c -' w (21) 
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where I(z) = 1 — 8z 2 + 8z 6 — z s — 2Az 4 £nz is a phase space factor (z=m c /mb) 
and (nib) is given by 

Cf f (m b ) = C 2 (M W ) + r,%C r (M w ) + _ V ^)C 8 (M W ) (22) 

Here Cg(Myy) is the co-efficient of the operator Q 8 = (g 3 /f 67r 2 )mf,S/j(T A " y T j4 6i 
G^ v in the renormalization group evolution, T A and G A are the gluon gen- 
erators and field strengths, and Ci represents the operator mixing with the 
4-quark operators in the renormalization group evolution— 

The Standard Model leading order prediction givefiJ BR[B — > X s ^f] = 
(2.9 ± 0.8) x 1(T 4 for m f = 174 GeV. Several analyses of b -> s + 7 in SUSY 
models within the leading order(LO)approximation have been carried out, and 
their effects on dark matter analysed. Generally one finds that the b — ► s + 7 
constraint is very stringent for fi > but less so for the case fi < 0. Thus 
for /i > one finds that as much as 2/3 of the parameter space', is. eliminated 
while for fi < only about 1/5 of the parameter is eliminatecE30SE3. More 
recently, the next to leading order corrections to the Standard Model iiave 
been computecO. Analyses including—these, give results which are similar to 
the ones based on the LO correctionfc3E2l'E3'Ej. 

5 Detection of Dark Matter 

A variety of possibilities for the detection of dark matter which involves more 
than just its gravitational effects have been discussed in the literature. These 
include :(i) scattering by nuclei of terrestial detectors ,(ii) annihilation in the 
centre of sun or earth, (iii) annihilation in the halos of galaxies, (iv) atomic 
excitations, and (v) superheated microbubbles. Presently, the most promis- 
ing of the above are (i) and (ii), and we shall discuss specifically (i) in detail 
here. Case(ii) is based on the idea that neutralinos in the solar neighborhood 
will be attracted to the center of the sun, where over a period of time they 
will lose energy and accumulate. Annihilation of these will produce neutrinos 
which in the terrestial detection apparatus will show up as upward moving 
muons(via their neutral current interaction). A similar situation is expected 
for the annihilation in the center of the earth. Background for these processes 
are significantly reduced _hecause of the small angular apertures around the 
center of sun and eartlcj'tj. Signals of this type can be picked up by neutrino 
telescopes. In case (iii) one expects that neutralinos in the halos of galaxies will 
annihilate into (j,e + ,p, ..)+X . However, the signal to background ratio does 
not appear very promising in this case. Case(iv) involves scattering of neu.- 
tralinos by electrons bound in atoms which may result in atomic excitationc3. 
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Realistic analyses of this possibility are currently lacking. Finally we men- 
tion the possiblity of superheated .xaicrobubbles as a cold dark matter detector 
which has been discussed recentljED. 

We turn now to a detailed discussion of neutralino-nucleus scattering^! C3. 
The fundamental interaction here is the neutralino-quark scattering. However, 
since the quarks are bound inside the nucleon and the nucleons are bound inside 
the nucleus, one has to deal with nuclear physics effects. We shall return to these 
later. We consider first the fundamental interaction involving neutralino and 
quarks. This interaction involves exchange of Z boson, squarks and Higgs. 
The dark matter in the Milky Way halo is estimated to have a velocity of 
v ~ 300fcm/s. In scattering involving these neutralinos the typical kinetic 
energy exchange is 0(lkeV M Xl /lGeV). Because the energy exchanges are 
small one can approximate the pole exchange diagrams by contact interactions. 
The total interaction is of the form 

L eff- L eff + L eff V 16 ) 

where L^fff is a spin dependent interaction and L^ 1 ) is a spin independent 
part. Lf^j} is given by 

L%f = {xil^Xi)[qiMLP L + A R P R )q] (24) 

where Pr.l are projection operators with Pr.l = (1 ± Js)/2 etc. Note that 
because of the Majorana nature of the neutralino xi7''Xi = and thus only 
the axial-vector bilinear involving the neutralino enters. The spin-dependent 
interaction receives contributions from the Z exchange and from the chirality 
diagonal (LL and RR) squark exchange terms. For the Z-exchange one finds 
that Af,Aft ex (n\ — n\) which means that this exchange not only depends 
on the Higgsino components of the neutralino [defined in Eq.(9)] but on their 
mismatch as well. On the other hand the chirality diagonal squark exchange 
depends on the gaugino components of the neutralino, i.e., one has 

^^K^ + Miy) 2 , A R ^{n 2 gJff (25) 

We note in passing that the spin dependent interaction discussed above in- 
volves no interference between the gaugino components and the Higgsino com- 
ponents. This aspect of the spin dependent interaction is in contrast to the 
scalar interaction which we discuss next. The effective scalar interaction in 
given by 

4?/ = (XiXi)(qCm g q) (26) 
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The scalar interaction receives contributions from Higgs exchange and from 
the squark exchange. The contribution to C from the Higgs exchange obeys 



C 



Higgs ^ 



gl 




I" COSOLH "I 

sinfi 




sinan | 
sm/3 I 


u — 


quark 




m h 


— sinaH 
I cos/3 J 




COSCtH | 

L cosfi ) _ 


d- 


quark 



(27) 



where h is the light CP even Higgs, H is the heavy CP even Higgs, Fh — (rii — 
n2tan6w){niCOsaH +n 3 sinaH), Fh = (ni —n,2tan9w) {n^sinan — n^cosan) 
and an is the rotation angle needed to diagonalize the h-H° mass matrix and 
is defined such thaO 



cos(2a H ) = =, tan2a H = - — „ 12 „ N 28 

where My is the Higgs mass matrix in the CP even sector. We emphasize here 
the fact that the inclusion of loop corrections t^thf Higgs masses is important 
in the computation of the Higgs mixing angl<M3'E!J. Jijurther, the contribution 
of the heavy Higgs exchange is in general not smaltHl. Thus while the heavy 
Higgs contribution is indeed small in the up quark sector, its contribution in 
the down quark sector can be very significant. The reason for this is that in 
the down quark sector the light Higgs contribution is proportional to sinan 
while the heavy Higgs contribution is proportional to cosan- Now over much 
of the parameter space one finds that tanan = 0(1/10). This leads to a 
suppression in the light Higgs sector and an enhancement in the heavy Higgs 
sector. Further, n^/n^ ~ tan(3 (see Appendix A) and hence Fh can become 
significantly larger than Fh. These enhancements in the heavy sector can 
sometimes even overcome the natural mass suppression of Tn\/m H making 
the heavy Higgs sector contribution quite large. Finally, the (chirality non- 
diagonal) squark exchange also contributes to the scalar interaction and one 
has: 



Qoc Tl TT^~T-2 n 3— -r-n-4 7T 29 

\ m q L ~ m xi ) l m |i ~ m xi ) S%Tl P COSp 

We note that the scalar interaction involves an interference of the gaugino 
and the Higgsino components of the neutralino. Thus if the lightest neu- 
tralino were to be strictly a Bino, the scalar interaction would vanish and the 
neutralino- nucleus scattering would be strictly governed by the spin depen- 
dent interaction. However, such a situation is never realised and the scalar 
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interaction always plays an important role in the event rate analysis. The 
scalar interaction is further enhanced by the fact that when one sums over all 
the quarks in the nucleus, one produces in the case of the scalar interaction a 
nuclear mass factor M^r, which enhances the contribution of the scalar term. 
Thus the scalar interaction is the dominant one for heavy nuclei. The event 
rate for a detector invoving a target nucleus of mass Afjv is then given by 



R = [Rsi + Rsd] 



Pxi 



O.SGeVcm- 3 



320fcm/s 



events 
kg da 



Here is the local mass density of Xlt and t^i is the incident \\ velocity. 
The spin dependent (SD) rate is given by 

Rsd = r 16W * lMjV l2 A 2 J(J + *) ' A ™ I" ( 31 ) 

where J is the nuclear spin and A is defined by < N | S ' i \ N >= A < N \ 
J | N >. For the spin independent part (SI) one has 

We note that for large Mn, Rsd ~ 1/Miv, while Ksi ~ Mn which implies that 
heavier the nucleus the more sensitive it is to the spin independent scattering. 

We discuss now the uncertainties associated with the computation of R$i 
and Rsd- 111 the computation of R$i most of the uncertainty is associated 
with the computation of the matrix elements of the operator m q qq between 
the nucleon states. One may write the matrix elements in the form 

< n\m q qq\n >= m n f { q n) (33) 

Estimates give /„ w 0.05, f s w 0.2 with as such much as 50% uncertainty 
in their determination^. Further the heavy quark contributions cannot be 
ignored. They couple to nucleons at the loop level via exchange of gluons and 
contribute to the mass of the nucleon. One may write < n\mQ QQ\n >—^m n 

(1- ^2 q fq 1 ^)- An uncertainty also exists in the determination of Rsd- This 
quantity depends on the matrix element 

< n | q^q | n >= 2s{ n) Aq (34) 
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Figure 1: Maximum and minimum of event rates for Ge with scalar interaction. 



Here s^ is the spin 4- vector of the nucleon, and Aq is the quark polarizability. 
Currently there are two determinations p£ Aq , one using the older EMC dataea, 
and the other using the new SMC dataa. The analyses using the EMC data 
and the SMC data give 



Am = 0.7 ± 0.08, Ad = -0.49 ± 0.08, As = -0.15 ± 0M(EMC) (35) 



Au = 0.83 ± 0.03, Ad = -0.43 ± 0.03, As = -1.0 ± Om(SMC) (36) 

Now the variations in Aq, and specifically the variation in As can give rise to a 
significant variation in Rsi- However, Rsi is typically not a large component 
of the total R except for light nuclei such as CaF2 . Thus the event rate analyses 
in lighter nuclei arc the ones most sensitive to variations in the values of Aq 
while the event rates in heavier nuclei are not very sensitive. 



6 Analysis of Event Rates for Universal Soft SUSY Breaking 

We discuss now the analysis of the event rates. Here we give the analysis under 
the assumption that the soft SUSY breaking parameters are universal at the 
GUT scale. As discussed in Sec 2, the parameter space of supergravity grand 
unification in this case is 4 dimensional and is spanned by 

mo,mg, A t , tan/3, sign(n) (37) 

As mentioned already in sec 2 the condition that ti not turn tachyonic places 
severe constraints on the parameter space. From eqs(12) and (13) one may see 
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Figure 2: Maximum and minimum of event rates for Ge with spin dependent interaction. 



this as follows: for A t negative A\ becomes large quickly since the A t term and 
the gluino term reinforce each other driving the t\ tachyonic. For A t positive 
one finds that A t also cannot get too large since that would make A\ large 
(rrig < ITeV) and again drive the light stop tachyonic. Numerical analysis 
shows that the no tachyonic condition limits A t to the following rangcE3E3 

- 1.0 ~ A t /m ~ 5.5 (38) 

Thus a heavy top mass eliminates a large amount of the SUSY parameter 
space. 

Event rates are computed under the constraint of the naturalness assump- 
tion on mo,irig < ITeV, and tan/3 < 25. A t is then restricted by the electro- 
weak symmetry breaking and by the condition that there be no tachyons,i.e., 
eq(38). In this range of the parameter space we compute the event rates for 
various target materials. The maximum and the minimum of the event rates 
for Ge are exhibited in Figs 1-3. Fig.l is for coherent scattering when we take 
into account only the scalar part of the interaction. Fig. 2 gives the analy- 
sis of incoherent part of the event rates where we take into account only the 
spin dependent part. Fig. 3 gives a composite of the event rates including the 
coherent part, the incoherent part and the total. A comparison of the maxi- 
mum and the minimum curves shows that for the case of Ge(a heavy nucleus) 
the maximum event rates are almost 100% dominated by the coherent part, 
while the minimum curves have a substantial component( as much as 40%) 
of incoherent part above a neutralino mass of about 60 GeV. Fig. 4 shows the 
corresponding curves for a light nucleus, CaF2. We also note that both the 
maximum and minimum event rates show a sharp dip in the vicinity of the 
Z-pole. This phenomenon arises as a consequence of the accurate method of 
computing the relic density which takes account of the correct thermal aver- 
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aging over the Z and the Higgs poles in the annihilation of the neutralinos(the 
lack of a second sharp dip in the event rates from the Higgs pole is due to the 
fact that one is integrating over the parameter space which smears the Higgs 
mass pole effects). The dependence of the event rates on the neutralino mass 
beyond the annihilation region involving the Z and the Higgs poles can also be 
understood. The event rates in this region are controlled by /i. Since [i 2 is an 
increasing function of -m^ , by Eq. (7) it is an increasing function of m-^ and 
m Q . As discussed above the maximum event rates are dominted by As/ which 
is proportional to (71-2713) [or (71277,4)] which by Eqs. (48-51) are the leading 
[0(Mz//i] terms. These decrease as \i increases. Thus we expect R to decrease 
as m-^-j and m Q increase. The analysis of Fig 3 bears this out. Finally we men- 
tion that the the maximum curve in the event rates has a strong dependence 
on tan/3. It increases as tan/3 increases. This can be seen from Eq.(27) where 
the leading d-quark contribution has C 1 Higgs ~1/ 1 cos(5^tan(5. 

7 Effects of Non-universal Soft Breaking 

Most of the previous analyses in supergravity grand unification have been 
done using universal boundary conditions for the soft SUSY breaking masses 
at the scale Mq. However, the framework of supergravity unification allows 
for the presence of non-universalities from a general Kahler potentialea. In the 
present analysis we wish to examine the effect of non-universalities on the pre- 
dictions on dark matter event rates. In the present analysis for specificity we 
shall rnpfipe -ourselves to consideration of non-universalities only in the Higgs 
scctoEjOca'Ea Thus we assume that the only non-universalities we have at 
the GUT scale Mq are given by 



where 5\ y 2 represent the deviation from the universal reference mass mo- All 
the remaining parameters are assumed universal. The important parameters 
that enter in the low energy domain are /z 2 , m. A , and m^ . The role of m^ has 
already been discussed in Sec. 2 We exhibit here the effects of non-universalities 
on the parameters /i 2 and m 2 A . We find 




ml(l + 5 x ); m 2 H2 



m 2 (l + S 2 ), 



(39) 




(40) 



Here 



Cx 



t 2 - 1 



1 



(1 
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2 




1 



1 



Do + 1 
2 



S 2 t 2 ) 



(41) 
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(deltal,delta2) = (0,0):tot(solid).coli(dots),iiic(da5hes) 
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Figure 3: A comparison of the maximum and the minimum of event rates for Ge for three 
cases: (i) the coherent part arising from the scalar interaction (dotted), (ii) the incoherent 
part arising from the spin dependent intcraction(dashcd), and (iii) sum of cases(i) and (ii) 
(solid). The maximum dotted line overlaps the maximum solid linc( sec Fig 1 for comparison 
of the dotted and solid lines) 



C2 = C 3 = -^(9 - t 2 e), C 4 = -±— x f (42) 

where t = tanft and where the functions e,f,g,k are as defined in ref(10). Sim- 
ilarly, for m 2 A one has 

m\ = ralD x + A\D 2 + m\D 3 + A miD 4 - \m\ (43) 

2 2 2 



where 



D i = — J C 1 g ' ~ 1 ( 1 5 — 2 ' " 



t2 + l i n t2 + l l \ n * 2 + l, (A ,, 

D 2 = -^jk, D 3 = --z—jig - e), D 4 = --2—j.f (45) 

In the above Dq defines the Landau pole in the top quark Yukawa coupling yo, 
i.e, 

y* n 1 a F W (ar\ 

Vo = mWo ; Do = 1 ~ 6yt W) ( 6) 
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(deltal,delta2) = (0 1 0);total(solid),coh(dots),inc(dashes) 
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Figure 4: A comparison of the maximum and minimum of event rates for CaF^ for three 
cases: (i) the coherent part arising from the scalar interaction (dotted), (ii) the incoherent 
part arising from the spin dependent interaction(dashed), and (iii) sum of cases(i) and (ii) 
(solid) 



and 

E(t) = (l + /? 3 )*(l + /? 2 )*(l + /?i)* (47) 

where Pi = ai(Q)bi/A-K : cti(0) are the gauge coupling constants at the GUT 
scale with ot\ = (5/3)ay, hi are the one loop beta function co-efficients defined 
by (6i,6 2 , 63) - (33/5, 1, -3), and F(t) = j* E(t)dt. 

Normally we expect |#i,2| < 0.5, and a similar level of non-universality 
in the relevant mass spectra at the electro-weak scale. However, under special 
circumstances enhancement of non-universalities can occur when the universal 
terms are suppressed. To exhibit this phenonmcnon wc consider the solution 
for /j, 2 . The dominant terms here are the first two terms, i.e., the terms which 
involve C\ and C2. Now C\ contains a universal part and a non- universal part. 
For large tan/3, the universal part accidentally cancels when Do=l/3 which oc- 
curs for a value of m t ~ 168 GcV, and consequently the non- universality effects 
get enhanced. An enhancement of non-universalities also occurs from the 
term, but in a somewhat different manner. Here the non-universalities are en- 
hanced when the residue of the Landau pole in A a vanishes. This happens 
when A t «0.6 rag. 
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Figure 5: Maximum and minimum of event rates for Xenon for the universal case when 
<5i = = 8 2 

8 Event Rates With Non-universal Soft SUSY Breaking 

In this section we analyse the event rates for Xenon with both universal as well 
as with non-universal soft SUSY breaking terms. Fig. 5 exhibits the maximum 
and the minimum of event rates for Xenon for the case 81 = = 82 . The 
general features of this case are similar to the case of Ge discussed earlier and 
the observations made there apply here as well. Next we discuss the non- 
universal case. When 82 > and Si < 0, the non-universal terms in Eq.(40) 
decrease /z 2 and can in some cases drive n 2 negative, eliminating such points 
from the parameter space. This effect is larger for small tan/3. Fig. 6 exhibits 
this case for 5 2 = -81 = 1 where the major effect is to significantly raise the 
minimum event rate (for lower values of ) since it is the smallest tan/3 that 
gives rise to the minimum values of R. Fig. 7 exhibits the opposite effect, i.e., a 
decrease in minimum value of R when 82 = — 1, #1 = 1 • The decrease in the 
minimum arises since is increased and R is a falling function of /x. Fig 8 
gives for comparison the composite of the three cases discussed above. We note 
that as m Xl increases and hence m g increases one finds that Ar increases and 
hence the Landau pole contribution becomes more dominant. In this region 
the effect of the non-universalities tends to get washed out. Further the effect 
of the Z-pole disappears for m Xl > 65 GeV. This explains in part why the 
minimum and the maximum event rates in the region of large gluino masses 
are insensitive to the effect of the non-universalities. Also, as m t is increased, 
Do becomes smaller and the Landau pole term in Eq. (40) dominates, making 
the effects of non-universal terms again less significant. Numerical analysis 
supports this conclusion. 

Finally we discuss the consequences of a more accurate determination of 
ilh 2 . One expects such a situation to occur in the next round of COBE like sat- 
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Figure 6: Maximum and minimum of event rates for Xenon for the non-universal case when 
<5i = -1 = -Si 



telite experiments which may determine the Hubble parameter to within 10% 
or better accuracy.This determination of the Hubble parameter will strongly 
reduce the error in flh 2 . Here we analyse what such an accurate determination 
will do to the prediction of the neutralino event rates. For specificity let us 
assume that tth 2 is constrained to lie in the narrow band 0.225 < Slh 2 < 0.275. 
The narrowing of the range of flh 2 more stringently limits the parameter space 
of the SUGRA model. The result of the event rate analysis for Xenon under 
the above constraint is shown in Fig. 9. An interesting new phenomenon is 
the appearance of peaks in the minimum curves. These sharp peaks originate 
from the constraints to prevent charge and color breaking minima£J(CCB con- 
straints). The position of the peak is determined by the mass parameter rnjj 2 . 
The peak moves to the right as m? H increases and moves to the left as m 2 H2 
decreases. The first situation occurs for the case <5i < 0, 8 2 > , and the second 
for the case 8\ > 0, 62 < 0. Numerical analysis for these cases as given by the 
dashed curve and the dotted curve support the analytical result. Comparison 
of Figs. 8 and 9 shows that a narrower range in Q,h 2 more sharply constrains 
the SUGRA parameter space and the allowed band of event rates shrinks. 

9 CONCLUDING REMARKS 

In this review we have given an analysis of event rates for the observation of 
dark matter in terrestrial dark matter detectors. Supergravity GUT models 
under the assumption of R parity invariance predict the neutralino to be cold 
dark matter over most of the parameter space of the model. The analysis 
of the event rates is carried out by ranging over the full parameter space of 
the SUGRA model within the prescribed naturalness constraints. Then we 
finds that the event rates can range from R = (10~ 5 to ~ 5) events/kg da. 
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Figure 7: Maximum and minimum of event rates for Xenon for the non-universal case when 
<5i = 1 = -S 2 



We find that for a heavy target such as Ge or Xenon the maximum event 
rates are almost totally coherent (which is controlled by the scalar interaction) 
while the minimum event rates can have a significant incoherent component 
(which is controlled by the spin dependent interaction) . For a light target such 
as CaF2, one finds that the minimum event rates are essentially incoherent 
while the maximum event rates have a very significant component of coherent 
scattering. 

We have also analysed in this paper the effects of non-universal soft SUSY 
breaking parameters in the Higgs sector on event rates. It is found that the 
non-universal soft breaking effects are most significant in the neutralino mass 
range ~ (35-65) GeV. In this region one finds that the minimum event 
rates can be decreased or increased by a factor of 10 depending on the nature 
of non-universality. 

The above analyses show that experiments with a sensitivity of R > 0.005 
events/kg da can explore a significant part of the parameter space. However, 
a further significant improvement in sensitivity by a factor of about 100 or 
more would be needed to explore the entire parameter space. For the case of 
non-universalities with 82 = —Si = 1 there is an improvement by a factor of 
about 10 in the event rates in the energy domain M% < 65 GeV. However, 
for the case 82 = — Si = — 1 one finds a reduction in the event rates by a factor 
of about 10 in the same mass region. Thus whether or not the event rates are 
enhanced or reduced depends on the nature of the non-universality. 

The above results have important implications for dark matter detec- 
tors. Thus the analysis shows that the parameter space of the model with 
62 = — Si = 1 can be exhausted completely in the mass domain < 65 

GeV by an improvement in the detector sensitivity by a factor of about 100. 
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Figure 8: Comparison of the maximum and minimum of event rates for Xenon for three 
cases: (i) universal case when 81 = = 62 (solid); (ii)(i)non- universal case when Si = — 1 = 
— 52(dashed); (iii) non-universal case when 5\ = 1 = — ^(dotted) 



However, much greater sensitivity will be needed to do the same for the case 
<5 2 = —Si = 1. Thus new detection techniques which have the potential of 
dramatically improving the detection sensitivities E3H, are important for a 
complete exploration of the parameter space of the supergravity unified model. 
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Appendix A 

To second order perturbation theory one finds El when (Mz//J>) 2 << 1: 



IM Z 
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M 7 



n-2 



= 1 



2 /i (1 — m^//i 2 ) rri2 — rh 
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- — sin29w 



sin2(3 - 
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2 fi 2 {l-m(/n 2 ) 2 



sin 8w 



-sin2j3 + 
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Figure 9: Comparison of the maximum and minimum of event rates for Xenon for three cases 
when 0.225 < (ih 2 < 0.275 :(i) universal case when 5i = = 62 (solid); (ii)non- universal case 
when <5i = —1 = — <52(dashed); (iii) non-universal case when &\ = 1 = — ^(dotted) 
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Eqs. (48-51) differ from the numerical computer results by amounts Srii ~ 0.03 
over a large part of the parameter space and are often significantly better. 



Appendix B 

As discussed in Sec 3, the conventional approximation in the computation of 
relic density has been the expansion of av in a power series expansion in v. 
However, in the annilation of the neutralinos one encounters Z and Higgs poles. 
A power series expansion of av in this case breaks down. One must therefore 
treat -the pole correctly in the thermal averaging and in the computation of 
J(xff3. We illustrate this below for the case of the Higgs pole. Here one has 



A 



av 



Higgs V 

((vi-e h y +1 l)) 



(52) 
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e h = (m 2 h -4m 2 Xl )/m 2 Xl (53) 
■y h = m h T h /m Xl (54) 

In the above m,h is the Higgs mass and Th is the Higgs width. J (xf) involves 
a double integration over a pole which for accuracy requires a high precision 
numerical computation. It is possible to significantly simplify the expression 
for J(xf) by carrying out one of the integrations analytically. One may then 
write J(xf) in the following form 

Jm ggs {xf) = f% wg ! {hh + ^rhh) (55) 



2V2m Xi y 7fc 



2 , „,2 



hH = \ dyy-le-viogl ^*';^ '* ] (56) 
z Jo e h + fh 



hH = \ dyy-h-y[tan-\ ^ ~ ^ + * ) + tan^)] (57) 

* JO ' lh lh 

The numerical analysis of J(xf) is significantly simplified using the reduced 
form. A similar treatment of the Z pole term can also be carried out. 
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